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Abstract 



The hadronic hght-by-hght scattering contribution to muon g — 2 is examined 



, based on the low energy effective theories of QCD, the Nambu-Jona-Lasinio 

, model and hidden local symmetry approach, supplemented by a general infor- 

mation concerning the asymptotic behavior of QCD. Our result is —52 x 10 

ix. 

D , with an uncertainty of ±18 x 10 which includes our best estimate of 

model dependence. This is within the expected measurement uncertainty 
of 40 X 10^^^ in the forthcoming experiment at Brookhaven National Labo- 
ratory. Our result removes one of the main theoretical obstacles in verifying 

the existence of the weak contribution to the muon g — 2. 
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I. INTRODUCTION 



A substantial improvement in the measurement of the muon anomalous magnetic moment 
= I [g^ — 2) is planned at the Brookhaven National Laboratory. The precision of the 
measurement is expected to reach the level of 

40 X 10"". (1.1) 

This is about 20 times more accurate than the best value available at present 

a^(exp) = 1 165 923 (8.5) x 10"^ (1.2) 

where the numerals in the parentheses represent the uncertainties in the last digits of the 
measured value. 

Compared with the electron anomaly, for which all contributions other than QED are 
negligible, the muon anomaly is more sensitive to shorter scales where the hadronic and 
weak interaction effects are important. Also, provided that the standard model prediction is 
known precisely, the muon anomaly will be a sensitive probe of physics beyond the standard 
model. A typical standard model prediction is 
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a^(th) = 116 591 877(176) x 10"". (1.3) 



We note that the uncertainty in ( |1.3| ) is comparable with the one-loop weak interaction 
correction 

a/,(weak-l) = 195 (1) x 10"". (1.4) 

Thus further improvement of the theoretical prediction is necessary in order to be able to 
confirm the existence of the weak correction term in a^. 



The uncertainty in ( |1.3| ) is dominated by the error associated with the estimate of the 
strong interaction correction to a^. The bulk of this effect is due to the hadronic vacuum 
polarization contribution, which starts at (9(q;^). (See Fig. 1 of Ref. |^ for the Feynman 
graphs which give this type of contribution.) Fortunately, this contribution is calculable 



without relying on our theoretical knowledge of strong interaction. The contribution 
to the a^(had.v.p.) can be expressed in the form 

amu\2 foo R(s)K(s) 



a^(had.v.p.) o{a^) = ^ — / „ ds (1.5) 



by applying the dispersion relation and the optical theorem. Here R{s) is the hadron pro- 
duction cross section in e^e^ collisions normalized to the lowest order formula for the fi'^fi' 
production cross section a^e'^e' — » /U^/x^) = 47ra^/3s. The formula (|1.5| ) enables us to 
reduce the issue of our ignorance of strong interaction dynamics to the experimental deter- 
mination of R{s) 0. The integral (|1.5|) has been evaluated by several groups For 



instance the estimate given in Ref. is 

a^(had.v.p.) = 7 068 (59) (164) x 10~", (1.6) 



where the first and second errors are statistical and systematic, respectively. Works PJTO 
which include more recent data are not too far off from ( |1.6| ), although the evaluation 
of uncertainties in the experimental data still varies considerably among authors. Future 
measurements at VEPP-2M, DA$NE,and BEPS are expected to reduce these uncertainties 
to the level of the upcoming experiment ( [I.ID JI(]|JTT| . 



On the other hand, the contribution of the hadronic light-by-light scattering diagram 
shown in Fig. |l| is potentially a source of more serious difficulty because it cannot be 
expressed in terms of experimentally accessible observables and hence must be evaluated 
by purely theoretical consideration. The purpose of this paper is to report on our attempt 
to estimate this hadronic light-by-light scattering contribution to the muon anomaly. A 
summary of our preliminary results has been given in Ref. [Q. We present the detailed 
analysis here. 

The paper is organized as follows. Sec. ^ starts with a survey of previously reported 
results on the hadronic light-by- light scattering contribution to the muon g — 2. With the 
help of chiral perturbation theory and Nambu-Jona-Lasinio (NJL) model, we find that the 
relevant diagrams associated with this contribution are the ones shown in Fig. |^ of Ref. 



13| . We also give an outline of strategies to solve the problems we have encountered. The 



next three sections are devoted to the treatment of the three types of diagrams; Sec. fT| to 
the charged pseudoscalar loop contribution of Fig. 0(a), Sec. ^ to the neutral pseudoscalar 
pole contribution of Fig. 0(b), and Sec. |V| to the quark loop contribution of Fig. 0(c). Sec. 
|VI| summarizes the present study and compares it with the recent result of Ref. [T^JT^ based 
on the extended Nambu-Jona-Lasinio (ENJL) model and discuss its implications. 



II. SURVEY AND IMPROVEMENTS 

This section begins with an overview of the previous studies on the hadronic light-by-light 
scattering contribution to the muon anomaly. We then point out a few problems associated 
with its evaluation, and describe the procedure which we have adopted to solve them. 



A. Previous Studies 

The muon anomalous magnetic moment receives important contributions from hadronic 
physics. Naive dimensional consideration suggests that the effect of the physics of the 
typical scale A higher than the muon mass is suppressed by (m^/A)^. This implies that 
contributions to from QCD will be dominated by nonperturbative aspects of QCD. Thus 
we are confronted with a calculational difficulty; the relevant hadronic contribution to the 
light-by-light scattering amplitude may not be calculable from first principles in the current 
stage of development of QCD. 

As the next best procedure, we may appeal to the chiral perturbation theory which 
attempts to describe the low energy dynamics of QCD in terms of hadrons. Its leading 
behavior is given unambiguously by the low energy theorems on the dynamics of pions (and 
kaons) which are the Nambu-Goldstone bosons resulting from spontaneous breakdown of 
chiral symmetry. The scalar QED calculation in Ref. [^] corresponds to the lowest-order 
evaluation in this context. Corrections to the lowest order results may be obtained by 
adding higher order terms of a power series expansion in momentum variables. 
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For the calculation of the muon g — 2, however, such a systematic chiral perturbation 
technique runs into some problems. Insertion of a vertex with high power of momentum 
into Feynman diagrams for the muon anomaly, as a correction to the hadronic light-by- 
light scattering amplitude, yields a divergent result. Thus we must resort to an alternative 
approach which unfortunately is more mo del- dependent. For instance, Ref. |^ introduced 
the vector meson resonances. It should be noted that the explicit incorporation of vector 
mesons allows one to compute higher order counterterms [l^ in the chiral Lagrangian. 



The resulting 0{p^) counter terms agree reasonably well with experimental determination. 
A well-known example of the success in this direction can be seen in the description of 
pion's electromagnetic form factor F{q'^), where q is the photon momentum. There, the 
vector meson dominance (VMD) model works even for q as large as the mass of p meson, 
Mp ~ 760 MeV . 

Now we shall return to our topic. From the point of view of chiral perturbation theory, 
pions will contribute to most significantly in the form of the diagrams shown in Fig. |^(a) 
and (b). A priori we do not know the magnitudes of photon momenta which are important 
for these contributions. For example one may attempt to estimate the contribution of Fig. 
0(a) in the lowest-order of chiral expansion which we will denote as a^(a, sQED). On the 
other hand we recognize that the VMD model describes the '7r+7r~7 coupling well for on-shell 
pions. Thus we are motivated to include the VMD model explicitly in the 71717 coupling. A 
naive approach, which leads to a^(a, nVMD), introduces vector meson to replace a photon 
propagator as ; 



I I Ml _ I I 



g2 q^Ml-q^ q^ - q^ 



(2.1) 



The numerical results obtained by following these procedures were |jT^ 

a^(a, sQED) = -0.043 7 (36) x (^^^ 

= -54.76 (46) X 10-", (2.2) 

and 



a.(a, nVMD) = -0.012 5 (19) x 

= -15.67 (2.38) X 10"", (2.3) 

respectively. We see a large modification when vector mesons are introduced. A natural 
question arising from this observation is : 

(Qi) Is the modification caused by the introduction of VMD model real ? If it is, why it 
seems to conflict with our expectation based on chiral perturbation theory that the vector 
meson effect is very small at low energies ? 

Next let us turn our attention to the diagram shown in Fig. |^(b). It includes the 7r°77 
vertex induced by the chiral anomaly. It is well-known that the effective interaction 

C = -^Tr^^^'^'F^^Fx^, (2.4) 

where fj^ ^ 93MeV is the pion decay constant and F^^, is the field strength of photon, 
describes the behavior of 71^77 vertex in the limit of zero pion momentum and on-shell 
photons. However, naive use of Eq. (|2.4]) for the 7r°77 vertices in the diagram of Fig. ^(b) 
leads to an ultra-violet divergent result. This is a signal that the interaction ( |2.4| ) is not 
applicable to photons and pions far off mass-shell and must be replaced there by some form 
factor. In Ref. |^ such a form factor was introduced by an ad hoc adoption of the VMD 
picture. Correcting a sign error in the previous calculation [^, this contribution was found 
to be 

a^{b) = -0.044 36 (2) x (^-^ 

= -55.60 (3) X 10"". (2.5) 

In the previous analysis the quark loop diagram in Fig. ^(c) has been treated as not 
independent of the first two diagrams. Rather it was used as an alternative approximation 
of the hadronic light-by-light scattering contribution to a^. If this assertion is correct, the 
result for quark loop p 

a^(c) = 0.048 (3) x (^^^ , 

= 62(3)xl0-", (2.6) 



in which constituent quark masses are used, should be nearly equal to the sum of (p73|) and 
( p.5| ). However, even their signs do not agree with each other. 

Therefore there arises the second question: 

{Q2) Are three diagrams shown in Fig. |^ independent after all ? 

We will examine these questions and explore the prescriptions for remedy in the next 
subsection. 



B. Improvements 

First consider the question Qi. It has been pointed out that the naive VMD model of 
IP] does not respect the Ward identities required from electromagnetic gauge symmetry [jl8 



We found further that it is not compatible with chiral symmetry. To solve these problems, 
it is useful to introduce VMD in a way that preserves chiral symmetry. This can be achieved 
by appealing to the hidden local symmetry (HLS) approach |T^. This formulation maintains 



gauge invariance and chiral symmetry explicitly and reproduces all the low energy theorems 
assured by chiral symmetry, such as the KSFR relation. The question Qi may thus be 
raised within the HLS framework. Keep in mind, however, that this approach is somewhat 
oversimplified. In particular it ignores higher resonances beyond the usual vector mesons. 
We must analyze and reevaluate the error in our final result taking account of the model- 
dependence. 

We shall now turn to the second question Q2. The previous work assumed that the 
quark loop calculation and the pion calculation are two distinct approximations to the same 
hadronic light-by- light scattering effect on a^. They should, therefore, yield the similar 



results and must not be added together. As was noted in Ref. |T^, however, in the extended 
Nambu-Jona-Lasinio (ENJL) model, the quark loop diagram contribution is independent of 
the other two so that all three contributions should be added altogether. 

This point can be made clearer by considering the 1/Nc expansion together with the 
chiral expansion. Table | lists the orders of each diagram shown in Fig. According to 
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the QCD-diagrammatic consideration, the pion loop needs at least two quark loops, and the 
pion pole diagram starts from a diagram in which at least one gluon propagates between 
a quark and an antiquark forming the pion. Thus a single quark loop contribution is not 
included in these two graphs. 

We may also examine this problem from the viewpoint of duality. In a dispersion relation 
for the light-by-light scattering amplitude, the integral over the quark loop diagram, from 
threshold to very high energies, is equal to the integral over the absorptive part due to all 
the hadronic intermediate states. Extension of this relation to local duality implies that 
the quark loop contribution approximates the hadronic contribution when certain averaging 
over a finite energy region is taken. Thus one may wonder if the quark loop diagram Fig. 
2(c), when embedded in the g — 2 graph, represents the entire hadronic contribution. 

We do not think so for the following reason. Consider a similar problem in the photon 
vacuum polarization diagram, and look at the cross section for low energy reaction e^e~ — >■ 
hadrons. Duality between the quark diagram and the resonances implies that 



where the integration range is sufficiently large so that the resonance is averaged out. In this 
sense, it is well known that the absorptive part of one quark loop diagram for the photon 
vacuum polarization is dual to the cross section from the vrvr threshold to about 1 GeV. This 
region is dominated by the p meson intermediate state, and the quark loop represents this 
contribution. Then there is a small continuum vrvr state contribution at lower energies. Near 
threshold, they can be calculated by chiral symmetry argument. Is this a part of the quark 
loop diagram or is it a non-resonating continuum? If it is a part of p, we suspect that a 
chiral invariant tttt interaction should be able to generate the p meson bound state. As is well 
known [^, however, the force between two vr mesons is not attractive enough to generate 
a bound p state. This supports the view that the vrvr intermediate state at low energies is 
independent of p resonance and hence independent of one-quark-loop contribution. For the 
photon vacuum polarization graph, pion loop graph and quark loop graph are independent 




(2.7) 
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and must be added together. 

The above argument suggests that the pion loop is independent of the quark loop in 
light-by- light scattering, too. The quark loop diagram corresponds to the sum of continuum 
hadronic channels as well as axial vector meson states. Of course, this is by no means a 
proof. Since it is impossible to prove this type of statements without solving QCD, we must 
keep this ambiguity in mind in our subsequent analysis. 

Let us summarize the above considerations and add a few corollaries: 

(1) The HLS approach avoids the inconsistency that has been observed in the naive VMD 
approach used in the previous analysis, i.e., violation of chiral symmetry and electromagnetic 
Ward identities (This will be demonstrated in Sec. |IIIB|) . 

(2) Three diagrams shown in Fig. || should be added. Especially the quark loop diagram 
which represents the averaged hadronic continuum effect in a certain energy region has been 
discussed as independent of the other two. 

(3) Contributions involving more loops of hadrons will be suppressed by a factor 
171^^/(4:71 f^,) compared to the contributions of Fig. |^(a) and Fig. 0(b). On the other hand, 
the 7] pole contribution from a diagrams similar to that of Fig. ^(b) may not be negligible. 
The magnitude of this contribution and the kaon loop contribution from diagrams similar 
to Fig. |^(a) deserves an explicit analysis. 

(4) As was mentioned already, naive use of Eq. (|2.4|) for the 71^77 vertices of Fig. ^(b) 
leads to ultra-violet divergence, indicating that ( p.4|) must be modified by a form factor 
far off mass-shell. Possible modification dictated by the asymptotic behavior of QCD will 
be discussed in Sec. VI. Here we simply note that the prescription adopted in Ref. in 
which the VMD was introduced merely as a convenient UV cut-off, can be justified within 
the HLS approach [^]. Note that the HLS Lagrangian can be obtained as an effective 
theory of the ENJL model [^]. Thus, a similar conclusion can also be reached in the 
ENJL model - the pion pole diagram contains two triangle loops of constituent quarks and 
p mesons is allowed to propagate before the quarks couple to photons. Fig. |^ shows this 
contribution diagrammatically. However, its evaluation needs some care, especially due to 
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the requirement of anomalous Ward identities |T^, as is described in Sec. [IV B . 



(5) In the quark loop diagram, the vector meson will affect the coupling of the quark to 
the photon. Using the ENJL model as a guide we determine the quark coupling to vector 
mesons. Its graphical expression is found in Fig. 4 of Ref. . 



III. CHARGED PSEUDOSCALAR LOOP 

A. Hidden Local Symmetry Approach 

For a complete description of HLS, the reader is referred to Ref. [|^. Mainly for the 
purpose of giving the Feynman rule relevant to our problem, we shall briefly discuss the 
formalism. The HLS incorporates vector mesons, such as p, as gauge particles of HLS, 
[SU(2)y]ioca\ in our case. The explicit form of the Lagrangian, assuming chiral symmetry 
[SU{2)l X 5'f/(2)/j]gj^j^^i and hidden local symmetry [SU{2)v]iocid is 

+ ^f^Bo {Tri^LM^n) + Tri^nM^^)} (3.1) 

where 

£^ = ^2Tr((af(x)f 
Cy = f^Tr {{a^Ax)?) . (3.2) 

In Eq. (|3.1| ) Fy and F^^, are the field strengths of the vector meson = fi'y^V^" (t'*, a = 
1,2,3, are Pauli matrices) and the photon A^, respectively, and gy represents the coupling 
constant associated with HLS. /^r is the pion decay constant (~ 93 MeV), and the coefficient 
a of Cy is an arbitrary constant to be fixed by experiment. consists of covariant 

derivatives of the basic objects and ^r(x) in the HLS approach: 



IIM 



2i 



dx, = (3.3) 



2i 
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where the covariant derivatives D^^l^r{x) are given by 



r3 



and ^/j contain the pion field 7r°(x) as well as the scalar triplet a°'{x): 

^^(x) = e^"(^)/^^e-*"(^)/^% (3.5) 
where n(x) = 7r"(x)^ and cr(x) = o""(x)^. The latter, on breaking the symmetry, will be 



absorbed into vector mesons V"^ to give them masses. In the last term of Eq. Bq is a 



dimension-one constant associated with the quark condensate [23 



Bo = -^{0\uu\0). (3.6) 

It combines with the current quark mass m„ in the mass matrix Ai = diag(m„, rrid) (we 

neglect the isospin violation due to the quark masses so that we set = henceforth) 
to give the pion masses 

ml± = mlo = 25om„. (3.7) 

In the unitary gauge a = for HLS, the relevant interaction terms for the present 
computation can be found as follows: 



+ (1 - a)eM^A^7r+7r- + 2e^p^,A^pj7r+7r-. (3.8) 
In this expression various masses and coupling constants are related to each other by |TP 



M'p = aglfl (3.9) 
9p = a9vfl, (3.10) 



gpirn = -agv, (3.11) 
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and Afj^ represents the photon field to the order e^. As is seen from Eq. ( p. 121) the complete 
vector meson dominance (namely g^T^.,^ = 0) is realized when a = 2. This is also close to the 
observed data. Note that Eq. ( |3.8| ) does not contain the p^p^-K^-K^ term. This is the crucial 
difference between the chiral Lagrangian ( |3.1| ) and the VMD model of Ref. 0. The absence 
of TT^TT^p^p^ coupling with no derivatives will be a common feature of chiral symmetric 
effective model, as implied by other models, too Il^j24| . 



B. Ward Identity 



Einhorn argued that the calculation of Ref. |Q does not satisfy the Ward identities 
among the couplings of vr and 7 required from the electromagnetic symmetry. The purpose 
of this subsection is to demonstrate its recovery in the present approach. For simplicity let 
us consider a 717 scattering amplitude and show explicitly that the relevant Ward identity is 
satisfied in the present approach. If we define the amputated Green functions G'^'^ and 
in momentum space by 



d%e'P^-y^ f d%e-'P^-y^ (0 |r [ji:^{z)A^{x)7c+{y,)n-iy,)] \ o) 



{2ny6\k-p^+p2)tT^{k;p,,p2) 



y2t 



T 



A'^{x)n+{y^)n-{y2) 



amp. 



amp. 



(3.13) 



and denote the (full) pion propagator as iD{p), the Ward identity can be written as 

- (l^Gf,^{q,k;pi,p2) = ^^'^^]^f' T^{k-pi + q,p2)- V^{k;pi,p2 - q). (3.14) 

In the naive VMD model of Ref. the photon propagator p5[ 



p 



(3.15) 



is replaced everywhere by 



—I 



Ml 



p2 p2_^2 p2[p2-Ml)' 



(3.16) 
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Thus, in this VMD model, the Green function (or G^u) defined above is simply obtained 
by the multiplication of one (or two) p propagator (s) to the corresponding quantity in scalar 
QED 



q2 _ Ml fc2 - Ml 



'2gnu - 



1 



{k +P2y - ml 



{2p2 + k)^{2pi + q)^ 



^^^{k■,Pl,P2) = e 



Ml 



{pi — ky — ml 

(pi +P2)m- 



(2pi - k),{2p2 - q), 



(3.17) 



Evidently the identity ( |3.14| ) cannot hold due to the difference in the numbers of p propa- 
gators between T ^ and G^,^. On the other hand, in the HLS approach, they are given, to 
the order of our interest, by 



G^,y{q,k;puP2) = -e 
1 



2l^g,,-^H^,ik)-^H^,iq)^ 
+ _ ^2 - ^H^pik)] {2p2 + kf [g,^ - ^H,^{q)] {2p, + qf 
_ _ ^2 [d'^P - \H.p(k)^ {2p, - kf [g,^ - \h,M)^ {2p2 - qT 



Tf,{k;pi,p2) =e[g^p- ^H^p{k) ) (^1+^2)'^- 



(3.18) 



where H^y{k) is defined by 



fc2 - M2 



^g^ivk k^ky^. 



(3.19) 



It is an easy algebraic task to confirm that the identity (|3.14|) holds now. In order to consider 
the recovery of Ward identity in more detail, we may add the k^k^/ M'^-ieim. to each p-meson 
propagator in Eq. ( p.l7|) by treating it as massive vector meson 

k^kiy 



M2 
p 



M2 
p 



g,^v 



M2 



gpv - Hf,„{k). 



(3.20) 



Then the expression ( p. 171) becomes 



G^u{q,k;p,,p2) = -e [2 {g^, - H^,{k) - H^,{q) + H,p{k)H'^^{q)} 

{gup - H,p{k)} {2p2 + kf {g^^ - H^^iq)} (2pi + qf 



{k + P2y - ml 
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{pi — ky — 

T,{k;p^,P2)=e{g,p~H^p{k)){p^+P2f. (3.21) 



The comparison of Eqs. (|3.21 ) and ( |3.18|) with a = 2 shows that the absence of the 



Hi3^{k)H^ Jq) term in ( 3.18| ) is responsible for the recovery of the identity (^Tl) . This 



is a consequence of the nonexistence of the direct p^p^n^-K -term, as has been stressed in 
Sec. |III A| . This argument apphes equally well to the light-by-light scattering amplitude 



caused by a charged pion loop. 

C. Muon Anomaly 

We can now evaluate contributions of diagrams of Fig. 2 to the muon anomaly a^. Let 
the vertex correction from a diagram 5* be denoted as A^(p, q) for the incoming photon 
momentum q, apart from the factor "ie" . Then the contribution to from the diagram S 
is given by 

a,{S)= hm Tr(P,(p,g)A^(p,g)), (3.22) 
where Pu{p, q) is the magnetic moment projection operator 

with the muon mass set equal to 1. The diagrams constructed from the interactions in 
Eq. (p.8|) can be classified in the similar manner as in Fig. 5 of Ref. 0. However, let us recall 
that the replacement ( p.l|) performed in the VMD model of works only if the p^p^TT^7r~ 
coupling term is present. In a theory with HLS ( restricted for simplicity to the case of 
complete vector meson dominance ( a = 2 ) ), however, there is no such term. This means 
that, while the replacement of the photon propagator ( |3.15| ) by ( |3.16| ) is performed as before 



if the photon line is connected to the pion through the 77r+7r^ coupling, the replacement 
must be carried out for either one of the lines but not for both, if the photon lines comes from 
a 777r+7r~ coupling. As a consequence, the contribution a^(HLS; A2) from the diagram in 
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Fig. ^ (which is topologically the same as the diagram A2 in Fig. 5 of Ref. |^), for instance, 
takes the form 



a^(HLS; A2) = a^(sQED; A2) - a^(sQED; {A^, 2)) - a^(sQED; {A^, 3)) - a^(sQED; {A2, 4)) 
+a^(sQED; (A2, {2, 3})) + a^(sQED; {A2, {3, 4})), (3.24) 

where a^(sQED; (A2, {2, 3, ■ ■ ■})) denotes the quantity obtained by replacing the photon 
propagators of the hnes 2, 3, ■ ■ ■ with the propagators of mass Mp. This differs from the 
calculation of Ref. by the absence of the terms 

+ a^(sQED; {A2, {2, 4})) - a^(sQED; {A2, {3, 2, 4})). (3.25) 

Since the 77?"'" vr" vertex receives no modification, the contributions of the diagrams Ci — C4 
in Ref. ||^ remain unaltered. 

The prescription for numerical evaluation of Feynman integrals follows that described in 
Ref. As in scalar QED, the B'^j which appears on the right-hand-side of Eq. (37) of Ref. 

must be changed to Bij. The correctness of this change can be shown explicitly in the 
same manner as in the Appendix B of Ref. . The renormalization |^ is required for cal- 



culating individual diagram since each diagram, not being gauge-invariant, has logarithmic 
divergence residing in the hadronic light-by- light scattering subdiagram. The evaluation of 



integrals is performed with the help of the Monte Carlo integration routine VEGAS [|28 |. 

Let us now summarize our results. To begin with we checked the scalar QED result in 
( p.2|) by writing new FORTRAN programs from the scratch. The result (obtained using 
VEGAS with 40 million sampling points per iteration and 60 iterations) is 

a^(a, sQED) = -0.035 57 (18) (^^^ 

= -44.58 (23) X 10"", (3.26) 

confirming the previous result in ( |2.2| ) but, of course, with much higher precision. The new 
evaluation of the p-meson contribution in the HLS approach yields (for 40 million sampling 
points per iteration and 60 iterations) 
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a^(a, HLS) = -0.003 55 (12) 



TT 



-4.45 (15) X 10"'\ (3.27) 



This result is about 3.5 times smaller than the VMD model result given in 

To see whether this reduction is real, we have evaluated the difference a^(a, HLS) — 
a^(a, nVMD) directly. The result is (for 40 million sampling points per iteration and 50 
iterations) 



a/,(a, HLS) - a^(a, nVMD) = 0.009 76 (4) ( - j 



12.23 (5) X (3.28) 



From (|3.27| ) and (|3.28|) we obtain 

a/,(a, nVMD) = -0.013 36 (14) 



= -16.74 (18) X 10"", (3.29) 



TT 
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which is consistent with ( [2.31 ). Thus the numerical works check out and the difference ( P3 
is real. Of course, the errors quoted above are those of numerical integration only and do 
not include estimates of model dependence. 



D. Discussion of Large Momentum Contribution 

As we have seen in the previous subsection, the p dominance structure has significant 
effects on the hadronic light-by-light scattering contribution to a^. In order to gain some in- 
sight in the dependence of a^{a\ HLS) on Mp and m^, let us introduce a function a^(a; m, M), 
where a^(a; HLS) = a^(a; m, M) for m = vtit^ and M = Mp, and examine it as a function of 
m and M numerically. Table ^ contains the result for m^r ( Mp ) dependence obtained by 
15 ( 30 ) iterations of Monte Carlo integration with one million sample points. From that 
table the following approximate asymptotic behaviors can be inferred: 
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TT 



a^(a; xm^, M) 4.75 x 10"^ x x^M - for x>3,M = oo, (3.30) 



a^{a; xm^, M) ~ +2.81 x 10"^ x x'^ for a; > 3, M = M^, (3.31) 



Tn,. \ f a^ '^ 



a^(a;m^,M) ~a^(a;m^,oo) + 0.23 (^^j (^-j for M > M^, (3.32) 

where a^(a; m^, oo) = a^(a, sQED). 

The results (|3.30|) and ( p.31| ) show that a^(a; m, M) depends on the loop mass m asymp- 
totically as for a wide range of "p mass" M. To appreciate the significance of these 
results, note that, for m » m^r, m in a^^a^m.,^^ M) — a^{a]m, M) may be regarded as a 
cut-off mass of the pion-loop momentum a la Pauli-Villars: The contribution of pion-loop 
momenta above m is suppressed in a^{a]mT^, M) — afj_{a;m, M). Thus the above depen- 
dences on result from the fact that the contribution of pion loop momenta larger than 
m drops off as m"^ as m increases. For instance, the contribution of pion-loop momentum 
higher than 800 MeV occupies only 7 percent of the total contribution (|3.26| ) or ( p.27| ). 



From these results, we speculate that the pion-loop light-by-light scattering amplitude 
(even with off shell photons) is governed by the region of small loop momenta carried by 
light hadrons. 

The result ( p.32| ) is inferred from the near constancy of M(a^(a; m,r; M) — a^{a; m.,^, oo)) 



for Mp < M < lOMp. This function decreases very slowly for larger M. Such an 
(instead of M"^) behavior seems to cast some doubt on the effectiveness of chiral pertur- 
bation theory since it implies that there is an appreciable contribution to 0^(0, HLS) from 
the region of photon momenta larger than Mp. To analyze this problem, let us recall that 
the customary argument in favor of the M"^ behavior is inferred from the fact that —M~^ 
is the dominant term in the p-meson propagator: 

1 1 p2 , , 

^ (3.33) 



for Ip^I <^ M^. As is readily seen by power counting, however, naive evaluation of the 
contribution of each term on the right-hand side of Eq. ( ^.331 ) to a^(a, HLS) leads to UV 



cut-off dependent results. In other words, the M ^ term has a divergent coefficient and 
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naive power counting argument fails. This is due to the fact that the Mp oo hmit and 
sub-integrations in the Feynman diagram do not commute. 

It is important to note, however, that the fact that photons of momenta larger than 
Mp contribute significantly to a^(a, HLS) does not necessarily mean the failure of chiral 
perturbation theory. The requirement on the photon momentum, which is a vector sum 
of pion momenta, can be considerably less strict insofar as the contribution to a^(a, HLS) 
comes mostly from small pion loop momenta. In fact this is what can be inferred from (p.30|) 
and (|3.31|) . For these reasons the behavior of ( p.32|) is not inconsistent with the chiral 
perturbation theory. 

For M = Mp, ( p.32|) can be written as 

ap{a, HLS) ~ ^-0.035 57 (18) + 0.23^^ (^^^^ , (3.34) 

where the first term is from ( p.26| ). a^(a, HLS) deviates from ( p.34| ) for Mp smaller than the 
physical value. It is seen from ( p.34| ) that the leading term happens to be nearly cancelled 
by the non-leading term for physical p mass. The smallness of the value ( ^.27| ) results from a 
(somewhat accidental) cancelation of ap{a, sQED) and the 0{mp_/ Mp) term for the physical 
value of the p meson mass. 

Before writing down our best estimate of the charged pion loop contribution to the muon 
g — 2, it must be recalled that the result ( p.27|) is based on the specific hadron model. In 



principle, any models which preserve chiral symmetry and the relevant Ward identities are 
the candidates for this computation. Any of these models is expected to lead to more or less 
the same hadronic light-by-light scattering amplitude provided that the chiral symmetry is 
intact. Eq. ( p.34| ) indicates, however, that the large photon momenta give rise to significant 
contribution to a^. Thus the hadronic structure in photon beyond the p mass may be 
non-negligible. The model dependence may enter here. 

In this computation, we have used the HLS approach. Even within this framework we 
assume further a complete p dominance. Also we could have chosen the version of HLS with 
higher resonances, such as Ai. All these would increase the uncertainty of our result. 



The previous work based on the ENJL model |jT3|] asserts that the QED result of pion loop 



should be included from the standpoint of systematic chiral expansion. As was mentioned 
previously, the ENJL model Lagrangian is always written in the form consistent with HLS. 
Then the pion loop contribution in that model reduces to the one obtained here when we 
approximate more complicated form factor of 7r+7r^77 and 7r+7r^7 PU|. 



In spite of these model dependences, we expect the total error to be within 20 % of 
the difference a^(a, HLS) — a^(a, sQED). This is because integrations over the photon and 
muon momenta are convergent in these diagrams and hence the contribution of large photon 
momenta does not distort our picture of low energy pion loop too severely. 

The kaon loop contribution is found to be about 4 % of the pion contribution. Taking 
these error estimates of the pion-loop and kaon-loop contributions into consideration, we 
present 



a^(a) = -0.003 6 (64) 



71 / 

= -4.5 (8.1) X 10-^\ (3.35) 

as our best estimate, including model dependence, for the contribution to of the charged 
pion loop part of the hadronic light-by-light scattering amplitude. Eq. (|3.35| ) replaces ( ^.3| ) 
obtained in Ref. [^] based on a naive vector meson dominance model which is inconsistent 
with the chiral symmetry. 

IV. NEUTRAL PSEUDOSCALAR POLE 



The first subsection here describes the detail of our prescription adopted in Ref. [|T2 
The second subsection describes extension of our method to include the pole-type axial 
contribution, and discusses the total contribution of this type. 
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A. Incorporation of Triangle Quark Loop 



For the purpose of later comparison, we record again the result (|2.5| ) for the neutral pion 
pole contribution a^(6) obtained in the HLS approach (for 5 million sampling points per 
iteration and 20 iterations): 

-55.60 (3) X 10"". (4.1) 



a^(6,HLS) = -0.044 36 (2) 



Here we used the newly written FORTRAN programs for evaluating this result. Note that 
a sign error in some part of the integrand in is corrected in ( |4.1| ). 

It is far from certain that the off-shell behavior, in particular, with respect to the pion 



momentum, is well-approximated by the use of the effective interaction (p.4|) modified by 
the HLS method. The examination of different off-shell extrapolation scheme will give some 
insight in the dependence of muon anomaly ci^(&) on the off-shell behavior. Here we choose 
the diagram shown in Fig. ^, which is again suggested by the ENJL model, as a model for 
such an extrapolation scheme, and evaluate it explicitly. 

Let us begin by noting that, in Fig. 3, the one-quark-loop subdiagram corresponding to 
7r°(g) — > 7(^1)7(^2) can be written as ^0 



TT/tt J Z2Z3P1 — ZSZ1P2 — ZiZ2q 

where [dz] = dzidz2dz36{l — {zi + Z2 + z^)). This amplitude is reduced to the one obtained 
from ( |2.4| ) in the limit pf = P2 = q"^ = ^, showing that it is normalized correctly. Note that 
the insertion of ( |4.2| ) into the Feynman diagram of yields a convergent result without 
recourse to the VMD model. If we use the notation a^{h]mT^, Mp^iriq) in analogy with 
a^{a;mT^,Mp) in the case of a^(a), such a contribution corresponds to a^(6; m^r, cxd, mg). 
The numerical evaluation can be carried out in a straightforward manner if we re-express 
( [4.2| ) in the form of momentum integral representation and use the general formalism |2^ 



for the evaluation of Feynman integral. The result is found to be (for 5 million sampling 
points per iteration and 20 iterations) 
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a„(6; m^, oo, m„) = —0.069 34 (5) (- 

VTT 

= -86.90 (7) X 10'^\ (4.3) 



where we have chosen the quark mass to be 300 MeV. The operator product expansion 
analysis on the short distance behavior of the amphtude A^is supports the use of constituent 
quark mass as niq in Eq. ( [4.2|) ||3T|] . 

The diagram in Fig. 0(b), with the VMD assumption added, can be calculated in a 
similar manner. As has been noted in Sec. [11 B|, the coupling of quark to vector meson is 



obtained with the help of the ENJL model. As a result we obtain (for 5 million sampling 
points per iteration and 20 iterations) 

3 



a 

(b) = -0.026 94(5) - 
Vtt 



= -33.76 (7) X 10-^\ (4.4) 

Here again, in order to examine what range of momentum governs a^{h), we perform the 
same analysis on a^{h;mT^^Mp,mq), where irtq is constituent quark mass of the triangular 
loop, as has been done for a^(a; m^r, Mp). Table |T| lists the results for the quoted quantity 
for various values of ( or Mp ) obtained by 20 ( or 10 ) iterations of integration with 5 
million sample points. For quark mass larger than 300 MeV, ap{b;mT^, Mp,mq) approaches 
a^(6; HLS) as m~'^, as is readily seen from the analytic expression for the triangle graph. For 
small rriq, it approaches zero. For instance, for rriq = 5 MeV, we have a^{b]mTj,Mp,mq) = 
—0.666 X 10~^(Q;/7r)^. The results in Table |T| are summarized by the following asymptotic 
form: 

ap{b; xm^, Mp, rUq) = -9.57 x 10"^ x x'^ (-^ for x > 3, (4.5) 



TT 



a^(6;m,,M,m,) = a^(6;m^,oo,m,)+0.3l(^j (^-) for M>3Mp. (4.6) 



These results show that the same consideration as in Sec. p.llD| also applies here. 



Note that the t] pole contribution, when the mixing among vr, t] and rj' is taken into 
account, amounts to 25 % of (^^ ): 
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af,{r] pole) = -0.005 29 (2) x 



= -7.305 (3) X 10^'\ (4.7) 

which is obtained using 5 milhon samphng points per iteration and 20 iterations. The p 
mass dependence of the rj contribution is hsted in Table ^ (for 5 million sampling points 
per iteration and 10 iterations). From this Table we obtain an approximate asymptotic 
formula 

a^(6;m^,M,mg) =a^(6;m^,oo, mg) + 0.11 (^) X (-) for M > 3Mp, (4.8) 



TT 



where 



a^(6; m^, oo, m^) = —0.020 08 (1) x 



TT 



5 



= -25.17 (2) X 10"'\ (4.9) 

This was obtained for 5 million sampling points per iteration and 20 iterations. 

Adding ( [4.4| ) and ( [4.7| ) and again estimating the model dependence to be within about 
20 % of the Mp-dependent term, we obtain 

a^{b) = -0.032 2 (66) (^^^ 

= -40.4 (8.3) X 10-^\ (4.10) 



This is the result for the pseudoscalar pole contribution given in Ref. [T^ • 

Further discussion of the off-mass-shell behavior of the 71^7*7* vertex is given in Sec. VI. 



B. Further Examination of the Pole Contribution 

After our summarizing paper, Ref. [|I2|, was submitted for publication, we learned that 
the hadronic light-by- light scattering contribution to the muon g — 2 has also been studied 
by another group in the large Nq limit within the framework of the ENJL model. Their 
initial result for the contribution corresponding to Fig. 2(b) disagreed strongly with our 
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result. Since then, however, it was found that this was due to a simple numerical oversight 



T^. Correction of this error brings their result closer to our a^{h). Nevertheless, their report 



stimulated our interest to study the axialvector pole contribution in some detail. 
In order to facilitate comparison with the results of |T^, we follow their method closely. 
In particular we use the notation which enables us to keep track of the normalization factor 
directly. For the operator js = qT^i'j^q, = qQ'^^q, where the isospin generator is 
normalized as 2tr (T^T'^) = 1, the PVV (pseudoscalar- vector- vector) three-point function is 
defined as 

njr(pi,P2) j d'x.e'P^-^^ J d'x^e'P^-^^ (0| Tj5(0)j?(xi)j,«(a;2) |0). (4.11) 

The part of n|^J^(pi,p2) corresponding to the 1-loop contribution, n^J^(pi,p2), is given by 
the well-known triangle loop graph 

2 1 

F{pl,pl,q') = l + h{pl,pl,q')-h{0,0,0), 

hipl,pl, q^) = 2m^ / dzidz2dz35{zi + Z2 + z^ - 1) ^ ' ~ ^\ 

J M^{Za) 
M'^iZa) =m^q- pIz2Z3 - pjz^Z^ - ^^2:1^2, (4.12) 

where is the momentum- cutoff which renders the quark loop contribution finite, ruq is 
the constituent quark mass, and r(n, x) is the incomplete gamma function 

POO 

T{n,x) = / dte-'e~\ (4.13) 
The leading I/Nq term of n|^J^(pi,p2) can be written as [32 



IV - ^g^,e,.„,p,P2^^^2(_^2) (^2(_^2) _ ^2) 

X ) - g^i-q') {l - Fiplpl,q')Liplpl)}], (4.14) 



where 



Lh' v') - mU-pI)mU-pI) .4 
"^^^-^^^ - {mU-pD -pIKmU-pD-pIY ^ ^ 
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and gs = Sn'^Gs/NcA'^ is the scalar coupling constant in the ENJL Lagrangian 



C 



(int.) 
ENJL 



NcM . 



X 1,3 



NcA^ 



(4.16) 



Here i and j represent the flavor indices, Nq is the number of colors and the parentheses as- 
sume the implicit sum over colors. The deflnition of various functions appearing in Eq.( [4.14[ ) 
can be found in Ref. [^. 

Next we turn our attention to the AVV (axialvector-vector-vector) three-point function 
defined by 

^f,uiPuP2) =^'l d'xre'^-^- j d^X2e'P^-^^ (0| Tj5.(0)jJ(xi)j,«(x2) |0) . (4.17) 

where js^ = qi^l^T'^q- A direct evaluation gives 

Ml(-g2) 



9A[~q 



t^z::^iPl,P2) 



-2mggA{-q^)j^^^-^^y—^tqaTliriPuP2) 
+2mggA{-q^) 2(_\\ _ 2 ^?c.n{^J^(Pi, P2) 



+2mggA{-q^ 



I 



m^[—q^) — q^ 



iqa n|^J^(pi,p2) 



+2mg{l - gA{-q^))- ^ 



^lu^{Pl,P2) 



Pi=0=pi=g' 



(4.18) 



where n^J^(pi,p2) is the 1-loop contribution of U^^^ {pi,p2)- Il^Api,p2) is a linearly 



divergent integral 



n^r>i,P2) = ^.^2/^(-l)tr 



7a75i-— 1 7/. 7 luj — u ' (^•^^) 

fr + f-^-niq fr-niq fr-f^-mq_ 

where Pauli-Villars regularization is understood. The last two terms of ( [4.19[ ) come from 

the presence of anomaly contribution when —iq^^^^{pi-iP2) is rewritten in terms of 

n^r(pi,P2): 



-^q^m:^{Pl.P2)L{plpl) = 2m, [ni^f^ (p.^pMp.^p^) - Ul^f^ip.^p^] 



p2=p2^g2=o 



(4.20) 
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a relation which was also used to derive Eq. ( [4.14| ). 

Now the contributions of pseudoscalar and axial vector intermediate states to the four- 
photon vertex graph can be written, for instance, as 



Ue 



+ 



-2i 



V 



NcAl 



AW a 



per 



(4.21) 



The pseudoscalar pole contribution can be extracted from ([4.21|) 



iAf,uiPi,P2] 



q2 _ ^2^_g2 



-iip^(p3,P4) 



where 



(4.22) 



Af,u{pi,P2) 



iiram 



" [l-gAi-q'){l-Fipl,pl,q')Lipl,pl)} 



(4.23) 



vr/.(-g2) 

For the following analysis, momentum dependences of various functions will be ignored: 
/2(-g2) ^ fl M^{-q^) ~ M2, (^^(-g') ^ 9a, - ^^(-g^) ~ A\q^ - ml) with A 
accounting for the wave function renormalization constant of pion 

dnil{-p^) 



A' 



Qp2 



(4.24) 



which is close ( and therefore set equal ) to unity, and is taken as oo. In this approximation 
Eq. (|]2|) reduces to 



iA^^{pi,P2)^ r-Zv4p^(p3,P4) 



(4.25) 



where the amplitude shown in (|4.2| ), multiplied by a function L(pi,p2) associated with vector 
meson dominance, is modified to 



A^,u{Pi,P2) = -^(^^^uPpPlP2FpyY{pl,pl,q^), 
Fpyy{plpl,q') = ^[l-g^{l- F{plplq')L{pl,pl))}. 



(4.26) 
(4.27) 
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The formal limit qa ^ 1 while keeping My{—q'^) at a fixed finite value reduces ( [4.26|) 



to ( |4.2| ) multiplied by L{p\,p'^. But such a limiting procedure is not self-consistent in the 
framework of the ENJL model. The term {l—gA) in ( [4.27] ) corresponds to the term necessary 



in order to recover the anomalous Ward identity missing in ( [4 .21 ) as claimed in Ref. P |. 

The contribution to the muon anomaly from the type of the graphs in Fig. |^ can be 
written as the magnetic moment projection ( see Sec. pi C| ) of 

le J s=i [2ti)'^ - f^- 

^ZJZJZJ ^2 ^ iAap{PuP2)iAxu{pz.q), (4.28) 

Pi P2 P3 Pi ~ ''''tt 

which includes the symmetry factor 2 and an approximation ~ 1. The internal lines are 
labeled according to Fig. The terms of Fpvv(Pi,P25 5'^) ( |4.27] ) consists of the point-like 



part, (1 —Qa), and the rest, gAF{p\,p\, q^)L{pl,p\). Thus the term proportional to (1 —QaY 
in the product of two A^^s in ( [4.2^ ) corresponds to the contribution which includes two 
point-like vertices which causes logarithmic divergence from the photon loop integration. 
This is handled by introducing the Feynman cutoff for the photon propagators 

q2 q2 q2_M^ q\M^ - q')' ^ ' ^ 

This procedure is formally the same as that used for incorporating the vector meson dom- 
inance property in (f4.1|) but with a new mass scale Mc instead of Mp. This allows us to 
check the program written for the present purpose. When Mc is set equal to Mp, the result 
( [4.1|) should be identical with the result of qa = 0, and the result (|4.4|) should correspond 
to that of qa = 1- This is explicitly confirmed by our program. We have also confirmed 
that the results corresponding to various values of qa { < qa ^ ^ ) always falls in the 
range between ( |4.4| ) and ( [4.1| ) for = Mp. In this way it is quite easy to observe that 
for any qa and M^, the cancelation among various terms cannot occur for the pseudoscalar 
pole contribution because all terms contribute with the same (negative) sign. Typical values 
of the TT^-pole contribution obtained for various values of qa and Mc are listed in Table |V| 
(10 million sampling points per iteration and 20 iterations for qa = 0.5, 2 million sampling 
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points per iteration and 20 iterations for the other). Although our calculation is not exactly 



identical with that of Ref. |jT^ since we have disregarded the momentum dependence of 
QAi—Q^)-, etc., our result should be approximately equal to theirs. It turned out that we 



were not able to reproduce the result in Ref. |T^ . 

The axial-vector meson contribution to four-photon vertex graph can also be extracted 
from (14211) 



{^^^^fKl:!{p^.P2)L{plpl)^ {^Ar:a'^nf;^{p„p,)L{plpl) 



—I 



^ ('4™t^{flr(P3.P-)i(Pl.rf) - nT(pM,,.,.,,j)- (4.30) 

From Eq. (3.44) of Ref. Qa/ Ia is found to be independent of Gy. Thus, all contributions 
in Eq. ( [4.30|) vanishes in the limit M4 —>■ 00 ( Gy ^ ), which, of course, should be the 
case. However the first term in Eq. (14.30|) may become numerically significant since its overall 
coefficient 

(4.3, 

where the equality is imposed by the ENJL model, is numerically large ~ 46 (for qa ~ 0.5 
IP^ ). Thus there remains a possibility that such a term contributes to the muon g — 2 with 
the same magnitude as pseudoscalar does, but with the opposite sign. 



An explicit calculation shows that (for gA = 0.5) the first term in ( [4.30|) , denoted as (1), 
and the second, denoted as (2), contribute respectively as 



a^{ai pole)[(l)] = -0.001 192 (1) x - 



TT 



a^(ai pole) [(2)] = -0.000 194 (2) x i^-j for = 1.0 GeV, (4.32) 

The first one was calculated by 8 million sampling points per iteration and 15 iterations, 
the second by 3 million sampling points per iteration and 15 iterations. 

We find that the axial-vector contribution has the same minus sign as the pseudoscalar 
one and is one order of magnitude smaller than the latter. Thus such a reduction of order 
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as was seen in Ref. cannot take place according to our calculation. In this respect we 



are in agreement with the corrected result in Ref. |]T5 . 

The axial-vector pole contribution ( [4.32| ) is negligible as a pole-type contribution, com- 
pared to the pion pole. Adding the new evaluation of the //-pole contribution (5 million 
sampling points per iteration and 20 iterations) 

af,{r] pole) = -0.011 69 (1) x (^^^ 

= -14.651 (5) X 10"" for = 0.5, = 1.0 GeV, (4.33) 

to the 77° pole contribution given in Table V for = 0.5 and Mc = 1.0 GeV, we obtain as 
the total pole contribution 



a^{b) = -0.045 9 (91) x i^-j 

= -57.5 (11.4) X 10"", (4.34) 

where the model dependence is again estimated to be within 20 % of the Mp-dependent 
term. This replaces the value in ( [4.10| ) as the total pole contribution. 



V. QUARK LOOP 

Inferred from the ENJL model, the quark loop diagram incorporating vector meson can 
be calculated by making the substitution ( p.l|) to photon propagators. This leads to 

a^(c) = 0.007 72 (31) (^^^ 

= 9.68 (39) X 10"". (5.1) 

To examine the quark mass dependence, we define afj_{c; xniq, M) with afj,{c;mq, Mp) = 
a^(c), where nig denotes the collection of such masses as m„ = = 300MeV and = 
500MeV and x the common scale factor. (Here we do not include c-quark contributions 
which have been included in the previous calculation without VMD. Note that the c-quark 
contribution in this case is negligibly small since the contribution of each quark of mass rriq is 

28 



then proportional to [|5[ . The contribution of c-quark in the present model will be found 
as further suppressed as is inferred from the mass dependence presented below.) Numerical 
studies similar to the previous ones are performed to examine quark mass dependence by 
iterating integration with one billion sampling points per iteration and 60 iterations. The 
pion mass dependence is also examined by iterating integration with one million sampling 
points per iteration and 50 iterations. 



The result is summarized in Table ^ and in the following asymptotic form: 



a^(c; XTTig, Mp) ~ 1.94 x 10"^ x x^^-" - for x > 3, (5.2) 



a^(c; mg,M) 



+0.044 - 0.43 



-4.0 ^ 


a 








)] 


\M 





3 



for M > 3Mp. (5.3) 



,7r/ 

Note that the suppression effect of vector meson is so large here that the value ( |5.1| ) is one 



order of magnitude smaller compared to ( p.6| ). However the strong damping property on 
the quark mass is consistent with the observation that only the physical degree of freedom 
is important at low energies [^. Algebraically such a rapid decrease occurs when all quark 
masses become comparable to Mp since the relevant mass scale of the system turns then to 
the quark masses so that the cancelation of the two terms in ( |2.1| ) begins. 

Again, we consider the errors arising from model-dependence to be within 20 % of the 
Mp dependent term. This is because integrations over the photon and muon momenta are 
convergent in these diagrams and hence the contribution of large photon momenta does not 
distort our picture of low energy quark loop too severely. We are thus led to 



ap{c) = 0.007 7 (88) 



TT 



9.7 (11.1) X 10-'\ (5.4) 



VI. SUMMARY AND DISCUSSION 



We have obtained the results ( p.35[ ), ( ^4.10| ) and ( |5.4| ) as the contributions of Fig. 0(a), 
(b) and (c), respectively. These diagrams have been discussed in Sec. |T| to contribute 
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most significantly and independently to the hadronic light-by-light scattering effect on muon 
anomaly a^, as guided by the use of chiral and l/N^. expansion. The Mp dependence of the 
contributions (|3.34| ), ( [4 .61 ) and ( p. 3D indicates that the integration over the photon momenta 



receives considerable contribution from the region where photons are far off shell. We have 
estimated that these high mass contributions should be well within 20% of the vector meson 
contribution, which leads to the large uncertainties assigned to ( p. 351) , ( [4.34|) and ( ^.41) . 
Combining these results we obtain 

(light-by-light) = -52 (18) x 10"^^ (6.1) 

This is almost within the error ( |1 . 1| ) in the upcoming experiment. Therefore, with the 
progress of measurement of R [|^, the accurate determination of muon anomaly by future 



experiment will actually show the presence of the weak interaction correction P,|3^ and 
serves as a new constraint on physics beyond the standard model. 

Let us now discuss possible causes of difference between our result and the recent result 
of Bijnens et al. []TB[, which is based on the ENJL model. For comparison's sake, let us list 



their results corresponding to Figs. 2(a), 2(b), and 2(c): 

a^,{a)BPP = (-14.5 22.8) x 10"^^ (6.2) 

for the cut-off fi ranging from 0.6 GeV to 4.0 GeV, and 

a^{b)BPP = (-72 186) x 10"^^ (6.3) 

and 

a^(c)BPP = (11.4 ~ 20.0) X 10"^^ (6.4) 
for the cut-off fi ranging from 0.7 GeV to 8.0 GeV. 



On the surface, the results of [|I^ seems to be more reliable than ours, being less depen- 
dent on assumptions outside of the ENJL model. On the other hand, their result is not free 
from ambiguities either mainly because their theory does not tell which cut-off should be 
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favored. In particular, it seems to be difficult to justify their results for large /i which lies 
well beyond the region of applicability of the ENJL model. 

The result ( p. 21 ) is about 4 times larger than our result (|3.35| ). This may partly be due to 



our simplifying assumptions, such as the complete vector meson dominance (a = 2) and the 
neglect of momentum dependence of various masses and effective coupling constants. Note, 
however, that the Lagrangian of [|l5l seems to have the 7r+7r^p°p° vertex. The presence of 
such a coupling (without derivatives) will be inconsistent with the low energy phenomenol- 
ogy. It also means that their Lagrangian does not satisfy the Ward identity (|3.14|) contrary 



to their assertion. In particular, their Lagrangian does not seem to incorporate the vector 
meson consistently, as is described in detail in Appendix A. If this is the case, it could 
explain the bulk of the difference. It should also be recalled that the smallness of our result 
( p.35| ) is a consequence of an accidental cancelation of two main terms for the physical p 



mass value. Such a delicate cancelation is not visible in the calculation of 15 



The contribution ( |6.3| ) is 2 to 5 times larger than our estimate ( |4.1CI| ). Since this is the 
largest term, it is the main source of disagreement between the two calculations. Actually, 
the low end value (—72 x 10^^^) of ( |6.3D is of the same order of magnitude as our value 
for a^(6; m^, oo, m^) given in ( [4.31 ). Recall that in the latter calculation the anomalous 
7r°7*7* vertex is approximated by a triangular loop of constituent quarks and photons are 
attached to the "bare" quark directly. If one assumes that any QCD modification softens this 
coupling, our result (^]3|) may be regarded as some sort of upper limit of the contribution of 
Fig. 2(b). On the other hand, the result (|6.3| ) increases with increasing cut-off /i beyond this 
"bound", suggesting that the result ( |6.3|) diverges logarithmically as p — > oo. This behavior 
is a consequence of the presence of a hard PVV vertex in their Lagrangian. Its prediction 
on the muon g — 2 must be viewed with severe reservation, however, since it is obtained 
by applying the ENJL model beyond its domain of validity determined by the cut-off A^^. 
In fact, such an unwarranted application of the model (with a hard anomaly term) violates 



unitarity as 7* goes far off shell and hence must be tempered with some form factor. In 
other words, any realistic theory must be consistent with unitarity, be it the ENJL model 
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or the HLS model. 

An examination of Fig. 3, in the hmit where both fermion triangles shrink to points, 
shows that the UV divergence arises from the integration domain in which the momenta 
carried by the photon 3 and pion 4 are small while the momenta carried by the photons 1 and 
2 are large. The far-off-shell structure of the 7r°7*7* vertex in such a region has been studied 
using the Bjorken- Johnson- Low theorem |3^, which shows behavior asymptotically, 



where g ~ gi ~ ^2 [l^- The case where only one of the photons (52) is far off-shell has also 
been studied by an operator product expansion technique Based on the latter analysis 
a formula of the form interpolating between Pi = and p\ = 00 

1 



l-(P?/(8vrV.^))' 

^ (6.5) 



1 - (p?/M2) ' 

has been suggested for the form factor F{p\,p\,q^) normalized similarly as that of ( [4.12[ ) 



The experimental data fit Eq. ( |6.5| ) very well with ~ (O.TTGeV/c)^ over the range 
2.0(GeV/c)2 to 20.0(GeV/c)2 In Ref. ^ it is argued that the off-shell behavior of 



the 7r°7*7* amplitude is represented reasonably well by the quark triangle amplitude ( [4. 21 ) if 
one takes account of the asymptotic freedom of QCD and a nonperturbative generation of 
constituent quark mass. The result of their analysis is consistent with those quoted above. 
These considerations suggest that our model based on Eq. ( [4.2|) may in fact be a reasonably 
good representation of the contribution of Fig. 2(b) ||40[] . 

There is relatively small difference between ( |6.4|) and ( p.4|) . The remaining difference is 
within the range of uncertainty caused by our simplifying assumptions. In fact the good 



agreement between ( |5.4| ) and ( |6.4| ) may even be an indication that we have overestimated 
the model dependence in ( p.4| ). As was mentioned already, this is consistent with the fact 
that integrations over the photon and muon momenta are convergent and do not distort low 
energy quark-loop picture too severely. 

It appears to be difficult to resolve the difference between our calculation and that of Ref. 



15| completely because of different approaches and because of the necessity to apply the low 
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energy effective theory of strong interaction beyond its safely tested domain. The complete 
resolution may have to wait for the lattice QCD calculation of the four-point function. With 
the rapid improvement of the computing power, such a day may not be too far off. 
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APPENDIX A: VECTOR MESON IN THE BPP CHIRAL LAGRANGIAN 



In this Appendix we show that the vector meson is not properly incorporated as a 



dynamical field in the chiral Lagrangian (5.6) of Ref. |jT5 



C = ^tr 
4 



p'^ V' 

9 



— tr [L^.L'^'' + R,,Rn . 



89' 



(Al) 



Here denotes vector meson field and A4 represents the quark mass matrix, Ai = 
diag(m„, m(i, m^), in the three-flavor case, and are the external vector and axial- 
vector fields respectively, and L^j^ and i?^,^ are given by 



■'flU 



d^igpu - dv) - d^igpf, - a^) + [{gp^ - a^), {gp^ - a^)] 



Rixv = d^{gp^ + ay) - dy{gp^ + a^) + [{gp^ + a^), {gp^ + a^)]. 



(A2) 
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The transformation properties of w^, a^, Pfj,, and the unitary matrix (U) consisting of the 
pseudoscalar meson are given by 

U^U' = VrUVI (A3) 

{v^ + a^) ^ {v'^ + a^) = Vr{v^ + a^)Vl + iVrS^VI (A4) 

K - «m) ^ K - = - + ^V^L^^Vi, (A5) 

{gpt. + a^) (^?Pm + O = + + ^V^R-^^^^, (A6) 

(^Pm - a^) ^ (ap^ - a'^,) = Viigp^ - a^)!^! + ^14 5^ ^2 (A7) 

for the chiral transformation {Vl, Vr). Then the covariant derivative P^f/: 

V^U = d^U - i(^p^ + a^)U + if/(^p^ - a^), (A8) 

in which vector meson p^ appears instead of to reahze the vector meson dominance, 
transforms covariantly. The vector meson mass term in ( |A1D is also chiral-invariant as 
shown below. (As a matter of fact it vanishes.) 

The transformation property of the combination {gpfj, — v^) is found to be 

(gp'f. - = (gp'f. + - K + «m) 

= VR{gp, - a,)Vl (A9) 



from (A4) and (|A6|). On the other hand, from (|A5|) and ( |A7|) , we find 



= VLigp^-a,)Vl (AlO) 

Since Vl and Vr are independent, we may consider the case where Vl = 1 and Vr is 
nontrivial. Then, from Eqs. (|A9D and (|A10|) we obtain 

(^Pm - v^) = VrIqp^ - Vf,)Vl (All) 

For simplicity, let us consider the two-flavor case, and set Vr = e*^'^^ (isospin rotation about 
the second axis). Then the RHS of eq. (|A1 1|) becomes 
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-{9P,-v,)r y*, (A12) 

where "*" denotes complex conjugation. 

Picking up the third component of both sides of ( |A11|) , for instance, we get 

{9P^-v,,f = -{gp,,-v^)\ (A13) 

that is, {gp^ - v^Y = 0. 

In a similar way we can prove that gp^ — v^ vanishes for other components. This means 
that is nothing but an external vector field and the vector meson has not been incorporated 
in the theory dynamical object. 
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TABLES 



TABLE L Orders with respect to l/N^ and chiral expansions of the diagrams shown in Fig. ^. 



Diagram 


l/Nc expansion 


Chiral expansion 


Fig. 1(a) 


1 


4 

P 


Fig. 1(b) 


Nc 


pS 


Fig. 1(c) 


Nc 




TABLE IL and Mp dependence of vr^ loop 
xa^{a; xttItt, M) for M = Mp and M = oo, and (xMp/rrip 


contribution. Table lists (xm-jv Im,,)'^ 
) X ra,,(a; rriTT, xM„) — a,,(a; TTiTT, oo)!. 


(vr/a)'^ X {xm-,^ /rripY 

X 

xap{a;xm.,r,Mp) 


(vr/a)^ X (xrH-Tr/m^)^ 
xap{a; xm^r, oo) 


X /7 . in' 771 T A/f ^ ) — n . { fl' 771 ^ CV) 1 


5 0.064 9 (44) 


-0.082 (8) 


0.234 (7) 


10 0.094 5 (45) 


-0.093 (9) 


0.190 (7) 


15 0.105 (4) 


-0.099 (14) 


0.155 (6) 


20 0.106 (5) 


-0.094 (17) 


0.141 (6) 


25 0.107 (5) 


-0.094 (18) 


0.124 (6) 


30 0.110 (5) 


-0.094 (25) 


0.120 (6) 
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TABLE III. m^r and Mp dependence of tt*^ pole contribution. Table lists (xm^/m^)^ 
xa,j,{b-xmT^,Mp,mq) and {xMp/nip) x[a^{b;mT^,xMp,mq) -ap{b;mT^,oo,mq)] . 



X 


(vr/a)^ • (a;m^/m^)2x 
xmT^,Mp,mq) 


/ / \3 f H IT 1 \ 

{-Klaf ■ {xMp/m^)x 
[an[b;mT^,xMp,mq) - a^(&;m^, (X),mg}J 


5 


-0.168 4 (1) 


0.336 7 (2) 


in 


—0 903 4 fTl 


976 ('?\ 


15 


-0.213 8 (1) 


0.235 9 (2) 


20 


-0.218 2 (2) 


0.207 3 (2) 


25 


-0.220 2 (2) 


0.186 7 (2) 


oU 


— U.zzi (^zj 


U.i (xj I \Z) 


lAoL/H/ iV. 

x[ai^[o,mT^,xMp, 


Mp dependence of 
'nT'q) —aij,{b;mTr,oo,mq)] . 


Tj pole contribution. Table lists (xMpjiriix) 


X 


(7r/a) 


• [xMp/mp,) X [ai^[b;mT^,xMp,mq) - a^(0; mTr, oo, m^jj 


5 




0.135 7 ((S) 


10 




0.114 5 (8) 


15 




0.108 6 (7) 


20 




0.087 3 (5) 


25 




0.078 9 (5) 


30 




0.072 (5) 
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TABLE V. 7r°-pole contribution for various values of gA and M^. All the values are listed with 
the factor (ck/tt)^ removed. 



9A 


Mc = 1 GeV 




= 2 GeV Me = 4 GeV 


0.3 


- 0.053 85 (2) 




- 0.089 58 (3) - 0.137 35 (4) 


0.5 


- 0.034 18 (1) 




- 0.057 06 (2) - 0.086 20 (2) 


0.8 


- 0.036 43 (3) 




- 0.044 03 (3) - 0.052 67 (3) 


V /I 

X tij 


j./vJDijii/ V i. ic^narK niciss aiiQ ivip 
u{c, xniq, ivip) ana [xivip/mfj,) x[a^(^c, 


dependence 
niq, XlVlp) — 


of quark loop contribution. Table lists 
an[c, mq,oo)\. 


X 


(7r/aj X X a^[c; xm 


q,Mp) 

xlO 


/ 1 \'\ furl \ 

(vr/a)'^ X [xMp/m^j) 
x[aij,{c;mq,xMp) - a^(c; m^, cx))J 





0.177 (6) 




— 0.489 (5) 


in 






— fl 48fi C^') 


15 


0.191 (51) 




- 0.460 (5) 


20 


0.193 (91) 




- 0.123 (1) 


25 


0.195 (142) 




- 0.402 (4) 


30 


0.194 (205) 




- 0.385 (4) 
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FIGURES 

FIG. 1. Hadronic light-by-light scattering (shown by the shaded blob) contribution to the 
muon anomaly. Solid line and dashed line represent muon and photon, respectively. 

FIG. 2. Representative diagrams which dominate the hadronic light-by-light effect on 
at low energies. Other diagrams are obtained by permutation of the photon legs, (a) Charged 
pseudoscalar diagram in which the dotted line corresponds to tt^, etc. (b) One of the tt*^ pole 
graphs, in which the dotted line corresponds to 7r° and the blob represents the 7r77 vertex, (c) 
Quark loop contribution, where quark is denoted by bold line. 

FIG. 3. Diagram of neutral pseudoscalar pole contribution with VMD and quark triangular 
loop. The bold dashed line represents the vector meson. The arrow attached to each internal line 
label indicates the direction of the corresponding momentum. Other diagrams are obtained by 
permutation of the photon legs. 

FIG. 4. A typical diagram contributing to aju(a). To facilitate correspondence with the text, 
a number is attached to each internal line. Other diagrams are obtained by permutation of the 
photon legs. 
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M. H. [There is no meaning in the above letters. It is inserted for printing out all pages 
in the two-sides. ( The total number of the files have to be odd, then. Sorry.] 
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FIG. 2 




FIG. 4 



